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discussed. A naive extrapolation is given of lattice results on the location of the critical end 
point to the continuum limit with the correct pion mass. An interpretation of current results 
on fluctuations is provided. A new method of analysis of experimental data is discussed, which 
provides a check on the removal of backgrounds, comparisons with QCD, and signals an approach 
to the critical end point. 
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1. Lattice 

When the light quark masses are finite, as they are in the real world, the critical point of QCD 
lies at finite temperature and chemical potential. As a result, lattice computations are beset with 
the sign problem. A method to bypass this and estimate the position of the critical end point (CEP) 
was given in ^ and used in [|2[ ||, ||]. The idea is to make a Taylor expansion of the pressure — 

TV 

where the Taylor coefficients, called the susceptibilities, are evaluated at \l% = 0, where there is 
no sign problem. The baryon number susceptibility is the second derivative of the pressure with 
respect to fig and has the expansion 

XB(T^ B )=^x {n+2) (T)^ B . (1.2) 

This susceptibility diverges at the CEP. The divergence can be diagnosed using the series coeffi- 
cients by the usual means and the CEP can be located. This is a physical point if the value of ju# at 
the end point is real. 

We implemented this on the lattice [Q, ||] using two flavours of light dynamical quarks. The 
quark mass was tuned such that the pion mass was 230 MeV. The lattice cutoff, A = 1 /a was varied 
between 800 MeV and 1200 MeV to estimate the range of lattice cutoff effects. All lattice compu- 
tations are done at finite spatial volumes. In this case the spatial box size was between 4 fm and 
6 fm near T c . As a result, the boxes were large in units of the pion Compton wavelength as well 
as the mean thermal wavelength. The temperature scale was set using three different renormaliza- 
tion schemes; the scale uncertainty is about 1%. The simulation algorithm was the R-algorithm. 
The main algorithmic parameter, the molecular dynamics time step was changed by one order of 
magnitude without any change in the results. 

Several issues remain to be addressed — 

1. What effect does the unquenching of the strange quark have? A claim that 2+1 flavour QCD 
has no CEP [§] turns out to be a lattice artifact, cured by decreasing the lattice spacing [g]. 
The numerical impact on the quark number susceptibility of unquenching the light quarks 
was earlier seen to be small |0J. This could indicate that unquenching the strange quark 
should not have a large effect on the position of the CEP. 

2. The state of the art is to use a pion mass of about 230 MeV. Decreasing this towards the 
physical value of 140 MeV should have a numerical impact on the prediction of the location 
of the CEP. A quantitative estimate can be made using results presented in [JsJ] . 

3. The global structure of the phase diagram may be more complicated. We have nothing to say 
about this. Current lattice computations address only the phase transition closest to Hb = 0. 

4. Finally, the series expansion is only carried out to a finite order (8-th order in our case). The 
limiting behaviour as the order is increased is intimately connected to finite volume effects. 
We discuss this next. 
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All lattice studies are necessarily performed at finite volume. The finite size scaling theory 
which is used to extrapolate to infinite volume is well established in the usual case where the 
simulation can be directly performed at the critical point. The maximum of the susceptibility then 
diverges as a (positive) power of the volume (a large effect), and the position of this maximum is 
shifted from its infinite volume limit with a (negative) power of the volume (a small effect). In this 
case, the simulation cannot be carried out at the critical point and these large and small effects have 
to be obtained (along with the position of the critical point) by analyzing the series coefficients. 

The result is simple. On any finite volume, V, the radius of convergence first seems to approach 
a finite limit, ^g(V), up to a finite order, Beyond this the radius of convergence will seem 

to diverge, since %b is finite at all finite V. The large effect is that becomes infinitely large 

as V — ► °°. The small effect is that jUg(V') changes by a small amount in the same limit. In the 
present day simulations the large effect is clearly visible, whereas the small effect is still hidden in 
the statistical uncertainties. 

Starting from present day lattice simulations, extrapolation to infinite volume, zero lattice 
spacing and the physical pion mass, all taken together would predict the most probable range 
for the location of the CEP to be T E = 165-175 MeV and jiif = 250-400 MeV. Note that there 
are many uncertainties and caveats in each of the extrapolations. An experimental search over a 
somewhat broader range of these parameters is therefore advisable. 



2. Experiment 

Away from a critical point the correlation length, £, , of baryon number fluctuations is finite. 
As a result, in any volume V there are ./V = V /£, 3 independently fluctuating sub- volumes. In the 
thermodynamic limit, as N — » °°, the fluctuations at temperature T are Gaussian — 

P(AB)oc eX pU^^-J, where 8B = B-(B). (2.1) 

One way to test whether the critical point is reached is to look for deviations from such Gaussian 
behaviour That such Gaussian behaviour should be observable is the content of JToj]. 

The current RHIC runs produce fireballs which freeze out in a region of the phase diagram 
which is not expected to contain the CEP. If lattice computations are correct about this, then one 
should see a clear signal of non-critical behaviour in present data. One way to analyze these is 
to construct the first few cumulants of the observed distribution, [B n ] for n < 4, and to extract 
from these the mean (B) = [B], the variance, a 2 = [B ], the skew, 5? = [B 3 ]/a 3 , and the Kurtosis 
= [B A ]/o A . At a normal point on the phase diagram one must have 

(B) oc V, oocy/v, y oc l /y/v, and X oc 1 /V. (2.2) 

In heavy-ion experiments the volume is not observable. So one is forced to use a proxy. The STAR 



experiment in a recent analysis Ql lfl used the number of participants as such a proxy; using this 
they verified the above power-law scalings. 

A clinching point would be to compare the microscopic cumulants with the QCD expectations: 

[B 2 ] = (T 3 V) ( ^ ) , [B 3 ] = (T 3 V) (^jrj > [B 4 ] = (T 3 V) X {4) ■ (2.3) 
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Figure 1: As T and fi are varied, the QCD predictions lie on a surface in the space of measurements 
(mi,m2,mj). The data will lie on this surface if all non-thermal behaviour has been properly removed. In 
the happiest case, this could yield a comparison with QCD and a measurement of T and pi. Deviations from 
the surface only in a small window of beam energy is a signal for the CEP. 



This is not possible until a few more questions are clarified. First, have all non-thermal sources of 
fluctuations (minijets, decays, etc.) been successfully removed? Answering this question needs a 
complete control over the systematics of the analysis. This has not yet been demonstrated. Next, 
at what stage of the evolution of the fireball were the fluctuations set up (what values of V and T 



should one use in eq. |2.1| )? This requires control over the theory of coupled hydrodynamics and 
diffusion [12]. 

In order to remove unmeasurable quantities like V and T from explicitly appearing in the 
measurements one has to construct different combinations of variables. Once all backgrounds and 
systematics are under control, the following combinations can be compared to QCD predictions — 

mi = yo 







1 


W] 


X (2) /T 2 


7-23' 


[*] 


x {4) 


2 


m 


- X (2) /T 2 - 


r 24 


[*] 


X [A) 


2 


[2?3] 
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m = -yO = 77^7 = .7^ ,„ = — • « 2 - 4 ) 

These three measurements can be used to extract ^( 2 > 3 - 4 ) under the assumption that all backgrounds 
have been removed and a comparison with lattice QCD predictions [Q] is possible. However this is 
a strong assumption. 

A neater analysis with less bias is possible. Plot the experimentally measured values of 
(m\,m2,m-s) in a three-dimensional plot. Lattice QCD predictions of these quantities can also 
be plotted in the same figure (see Figure [[]). As one varies T and /j.b, the theoretical predictions 
trace out a surface. If the experimental points do not lie on it, then non-thermal sources have not 
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been completely removed. On the other hand, if they do, then one can estimate the values of T and 
He the experiments correspond to by comparison with the QCD predictions. 

Near a critical point, one can never satisfy V/i§ 3> 1. As a result the central limit theorem will 
break down. In a static system one would expect the Kurtosis to diverge. However the expansion of 
the fireball rounds off the transition and renders B, finite [|l3|]. Nevertheless, E, has a maximum near 
the CEP. Since -J€ ~ t, 3 [jl4j] the Kurtosis peaks near the CEP. For the measurements suggested 
above, one would have 

mi - ^ 7 -^)/ 2 , m 2 oc^ 5 -^ m 3 - ^ 5 - r ))/ 2 . (2.5) 

In other words, all these quantities would have non-monotonic behaviour as a function of the beam 
energy if the beam-energy scan passes through the vicinity of the CEP. If desired, such analyses 
can be easily extended to higher cumulants (when the system deviates from a Gaussian, the higher 
cumulants become easier to measure). 

One can use the plot of Figure |l| in the search for the CEP. Tune the background subtraction 
and cuts so that the present data for (mi,m2,mj) lie on the QCD surface, as it should. Then in a 
beam energy scan, the data will lie on the surface whenever the system lies away from the CEP. 



In the vicinity of the CEP, however, effects such as those discussed in Q13| ] drive the system away 
from thermodynamic equilibrium. As a result, in this small window of energies the experimental 
data will deviate from the surface, signaling critical slowing down as a direct probe of the nearness 
of the CEP. 
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